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I. INTRODUCTION 

Today, the general theory of relativity is a generally recognized metric theory of gravity and underlies our knowledge 
of the structure of spacetime. This theory is based on Einsteins equivalence principle (EEP). One of the consequences 
of this principle is the gravitational shift of spectral lines (the so-called gravitational redshift), which is equivalent to 
the dependence of the clock rate on a localgravitational potential ip . According to the local positional invariance 
(LPI) principle, which is a part of the EEP [l|, the gravitational shift of the clock frequency is universal and must not 
depend on the type of clock used. On the other hand, attempts to develop a quantum theory of gravity sometimes 
suggest that the EEP should be violated at a certain level; thus, these attempts stimulate experimental verification 
of the EEP and, in particular, the LPI to an increasing degree of accuracy. The LPI principle can be tested by 
experiments on measuring the gravitational redshift. According to this effect, the dependence of the frequency v of 
any clock on the potential . in weak gravitational fields (ip / 1, c being the speed of light)is given by i/q (1 

+ if / c^), where vq is the frequency of a clock in the absence of the gravitational field (the proper frequency). If we 
assume that the LPI principle is violated, then this dependence can be represented as [l| 

l^A = '^o[l + {l + f3A)v/c'], (1) 

where the dimensionless quantity /3a characterizes the deviation from the redshift law, which follows from the EEP, 
and the subscript A indicates that the frequency va may depend on the specific type of clock used. 

There exist two methods for verifying the LPI principle. The first method (the classical experiment on the grav- 
itational redshift) is based on the use of identical clocks situated at points with different values of the gravitational 
potential ( ipi and (f2 ) and on the measurement of the frequency difference 

I'Al - VA2 = 1^0 [1 + (1 + Pa) (<^i - ^2) /c^] ■ (2) 

In f^l, the authors carried out numerous experiments to determine possible values of /3a with increasing accuracy. 
Comparing the oscillation frequencies of two separated hydrogen masers (H), they obtained the minimal upper limit 
\/3h\ < 7x 10^"'^ . The second method for verifying the LPI principle (a null gravitational redshift experiment [l|) 
is based on the comparison of the rates of nonidentical clocks A and B situated at points with identical values of 
the gravitational potential under time variations of the potential f (t) itself. In this case, the violation of the LPI 
principle leads to a nonzero frequency difference 



VA-i^B = i^oi/3A- l3B)v{t) /c^- (3) 

A large number of experimental studies devoted to the verification of the LPI principle by formula ©were based on 
the comparison of the frequencies of atomic transitions in different substances under the variation of the gravitational 
potential (p (t) due to the orbital motion of the Earth around the Sun. The least upper bound for the possible values 
of Pa — Pb was obtained by this method in [3|] by comparing the frequency variations of a hydrogen maser and 
the cesium frequency standard (Cs): \Ph — Pcs\ < 2.1x10"^ . In a null gravitational redshift experiment 

was carried out by comparing the frequency of an atomic transition with the eigenfrequency of an electromagnetic 
resonator; i.e., the frequencies to be compared were of completely different physical natures. In the measurement 
of the frequency variations of a hydrogen maser and of a superconducting electromagnetic resonator yielded an 
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upper bound of IPh — Pres\ < 1.7x10 ^ , whereas, in the comparison of the variations of the electron transition 
frequency in iodine molecules (I2 ) and the eigenfrequency of a cryogenic optical resonator yielded an upper bound of 

1/3/, -/3res| <4X 10-2 . 

Experiments with clocks of different physical natures yielded an upper bound that is three orders of magnitude 
greater than that obtained when verifying the LPI principle by comparing only atomic transition frequencies in 
different substances; this is associated with the large intrinsic noises of the experimental setups used in However, 
it would be desirable to obtain a maximally accurate estimate for (3a - Pb precisely in experiments with clocks of 
different physical natures. Below, we will show that this may provide an experimental verification of the relativistic 
generalization [6|-[3| of the classical theory of elasticity. 

In the present paper, we consider the possibility of verifying the LPI principle by a stationary horizontal double 
laser system Q one of resonators of which contains a cell with a nonlinear ly absorbing gas. 



II. OSCILLATION FREQUENCY OF A GAS LASER WITH A NONLINEARLY ABSORBING CELL IN 

A GRAVITATIONAL FIELD 

In [§] , the effect of the gravitational field of the Earth on the oscillation frequency of a gas laser was investigated on 
the basis of a covariant generalization of the Lamb theory. The investigation was carried out within a semiclassical 
approximation, when the electromagnetic and gravitational fields are described classically, while the active medium 
is described quantummechanically. The formalism of the Lamb theory [To| was applied to derive equations for a 
self-consistent field. As shown in [9|, on the one hand, the formation of a macroscopic electromagnetic field in the 
active medium is determined by Maxwell's equations, which are essentially relativistic and contain the gravitational 
field via metric coefficients and their derivatives; on the other hand, it is the mechanical subsystem, the resonator 
evolving in the gravitational field by different laws, that performs the selection of harmonics. 

In Q , oscillation equations for a single-mode linear laser situated horizontally on the Earth's surface were obtained 
with regard to Newton's gravitational potential (/? at a given point: 

^ + ^ - ^m] = -2^^ - ^) Cm, (4) 



dt 



dt 2 \ J V c2 

Here, E and are the slowly varying amplitude and phase of the generated electromagnetic field; c is the velocity of 
light; bj is the oscillation frequency; is the natural longitudinal eigenfrequency of the m the mode of the resonator 
in the gravitational field, 

= f^mO ( 1 ~ j ' ^rnO = ^mC = — ^ ; (6) 

Aflji is the linewidth of the empty resonator; and L ~Lq (1 +^</? / ) is the length of the resonator, where the 
phenomenological parameter ^ is determined by the solution of the elastodynamic problem of the evolution of the 
mechanical length of the resonator in the gravitational field and Lq is the resonator length in the absence of the 
gravitational field. In the presence of a nonlinearly absorbing gas in the resonator, the polarization coefficients 
and Sm can be represented as fTo'] Cm = (x'a + Xb)^ Sm— [Xa + Xb)^^ where the real ( x'a) ^"^^ imaginary [Xa ) 
parts of the nonlinear susceptibility of the active medium, calculated by the method of Q , are given by the following 
expressions up to the second- order terms with respect to the field (under the assumption that the homogeneous 
amplification linewidth is small compared with the Doppler width and the detuning of the oscillation frequency from 
the center of the amplification line is small): 

x; = da^^(i + 2c^) 

(7) 

X {l - [1 + (1 + ^] (Aa) CaE^] , 

Xa" = -da (l + {l-[l + D (A J] aaE^] , (8) 
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'*^™0 — Lo ■ 

Here, T)i2a is the matrix element corresponding to the electric dipole transition in the atoms of the active medium 
between operating levels 1 and 2, which are characterized (in the absence of the gravitational field) by the decay 
constants and 72^, respectively, and by the homogeneous linewidth F^; is the proper (independent of gravitation) 
frequency of the atomic transition [7]; Na is the population inversion caused by pumping (Na > 0); ks is the 
Boltzmann constant; and Ta and m^ are the temperature and the mass of the atoms of the active medium. 

Expressions for the real {x'b ) ^^'^ imaginary {xb" ) parts of the nonlinear susceptibility of the absorbing medium 
are obtained from formulas ([7]) and ([8]) by changing the subscript a to 6 (in this case, iVf, < 0). 

For a steady oscillation mode, we derive from ^ an equation for the intensity of the electric field. Using this 
equation and introducing a dimensionless quantity r = db/da > 0, we obtain from ([4]) the following equation for the 
oscillation frequency: 

UJ-n,n^ {uJa - UJ) {pa [l + (1 + f ] " QaD (A,)} 
- {tOa - Lo) {Pb [1 + (1 + e) f ] - qbD (Afc)} , 



P'^ ~ (l-r)(l-2aB^) ' 



la — 2r2(l-r)(l-2a£;2) ) 



(10) 



Pb ~ (l-r)(l-2a£;^) ' 

(11) 

2r<^{l-r){l-2aE^) ' 



a„[l + D(A„)]-ra,[l+J(Ai,)] 
2(l-r) 



(12) 



For a stable steady oscillation mode, it is required that a > 0. Since the atoms of the absorbing cell are usually 
characterized by a narrow absorption linewidth (F^ ^ Ta), we have g;,> qa', i.e., the main nonlinear phenomenon is 
the nonlinear pulling of the oscillation frequency to the center of the absorption line (frequency self-stabilization) (llj 
Indeed, setting A^.b ^ 1, we derive the following expression for the oscillation frequency from Eq. (|lip : 



u -ujb+ 1+s + -T+s- 



X {[l + (l+0 ^]Pa-qa}, 



(13) 



The quantity 



S 



i + (i + ^ 



{Va -Pb) - (qa - qb) 



(14) 
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is called a self-stabilization factor. For typical parameters of the absorbing cells (g^ ^ Pa,b)i the inequality S ^ qb 
1 holds; therefore, the oscillation frequency is determined by the absorption frequency lui, of the atoms of the cell in 
the gravitational field, 



'^t (l + J) • (15) 



In the case of a laser without an absorbing cell in the resonator, the oscillation frequency is determined (see formula 
(fTTjl for r = and Sa -^il) [9] by the eigenfrequency fimo of the resonator: 

uj = n^Q ( 1 + ^ ) • (16) 



III. COVARIANT PROPAGATION EQUATIONS FOR WAVES IN AN ELASTIC MEDIUM 

The eigenfrequency 

of a resonator depends on the geometrical dimension L(t) of the resonator; to determine this dimension in the 
gravitational field with potential ip(t), one should solve the elastodynamic problem, which predicts the value of the 
parameter ^. 

At the classical level within Newton's limit (i.e., in three-dimensional Euclidean space when the speed of light 
c— > oo), the equation for the propagation of elastic waves is well known. Within this theory, the displacement Ui of 
an element of a body under a volume force is described by the equation 

U. + F.= '-^. (18) 

Here and below, we assume that summation is performed over repeated indices; the Roman indices run over the 
values 1, 2, and 3, and the Greek indices run over 0, 1, 2, and 3; the dots denote differentiation with respect to time; 
aik is the stress tensor; and p is the material density of the medium. In the case of a gravitational filed, the force 
Fi can be represented as Fi = Viip. Hooke's law, which relates the stress tensor to the tensor of small deformations 
of the body £mn, is expressed as aik =C ikmnSmn in the classical theory; here, Cikmn are the elastic constants of the 
body and, due to Saint Venants condition for the consistency of deformations, the tensor of small deformations .mn 
is related to the vector of small displacements Ui as eij= {Uij + Uj^i)/2. Henceforth, Uij — dUi jdxK 

The classical theory of elasticity has a long history, whereas the development of the relativistic theory of elasticity 
began rather recently. The first attempt was made by Weber in the late 1950s [12 L The equations derived by him 
are used for describing the response of an elastic pig to a weak gravitational wave 13 |,[]3 . In its modern form, the 
relativistic theory of elasticity was constructed only in 1972 by Carter and Quintana jl5l |. This theory is still being 
developed (see, for example, [lB|j[I3|)- ^ significant contribution to the development of this field of science was made 
by Maugin. For instance, he proposed in [18j a relativistic equation for the propagation of elastic waves. 

To calculate the evolution of the resonator length L(t) in a variable gravitational field (see formula (HZ])), we 
apply the propagation equation for elastic waves [13|- Within this approach, Hooke's law and the equation of elastic 
vibrations in spacetime V ^ are expressed as 

ih — CikranSmn^ (1^) 

e'ifc = -<7{i\Li\k) + c^RiOko, i,k = 1, 2, 3. (20) 

The boundary conditions are written in the standard form, {(JikN^)]^: = Fi, where Fi is an external surface force 
of nongravitational origin and A'^*' is a normal vector to the surface of the body S. Everywhere below, we assume that 
forces of nongravitational origin do not act on the body, i.e., Fi — 0. 
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For a homogeneous isotropic medium, the elasticity tensor is represented as 

Ciknin — ^^ik^mn {^ini^kn ^in^km] ; (^-^) 

where X.fi, — const are the Lam'e coefficients related to the longitudinal and transversal velocities of elastic waves in 
the medium by the formulas 




a; = \ , at (22) 



Within this theory, in contrast to the classical one, the deformation of a medium consists of a background defor- 
mation hik/2 and a true deformation Uik-£ik = Uik + hik/2, where = Uij = {Ui,j + Uj j)/2. 

Consider the metric of a tidal action, i.e., the metric that describes a quasistatic uniform isotropic Newtonian 



gravitational field: 



dS^ ^ + c^dt^ + S.jdx'dx^. (23) 



The choice of the metric in this form has been motivated by the following facts. This metric has the form of the 
Schwarzschild metric expressed in isotropic coordinates up to terms on the order of M/r, where M is the mass of the 
gravitating center and r is the distance to this center. In this case, the function ip has the sense of the potential of the 
gravitational field at the observation point. Since the Earth moves along an elliptic orbit around the Sun, r = r(t) for 
a terrestrial observer; therefore, Lp = f{t); however, at any moment in time t, metric (I23p is the Schwarzschild metric 
expressed with the above-indicated accuracy. In addition, the variations of the gravitational potential at a given point 
of the Earth are also associated with the motions of Moon around the Earth. It is the varying gravitational potential 
at the observation point that we model by the function f{t). The values of this function at a given point of space 
at any moment can be produced to a high degree of accuracy. Note that, compared with the rates of processes in a 
laser, the variation rate of the potential . is extremely low. Therefore, for such processes, we can assume to a high 
degree of accuracy that the potential (p is constant. Further, we distinguish the constant {^o) and variable {Sip{t)) 
parts of this function: (p{t) — (po + 6ip{t). 

In this case, the elastodynamic equations (I20p can be rewritten as 

C/,fc = i [3A(5(,|,0,i|fc) + 2fiU^m^nk)] , (24) 

where Q = Uk,k- Equation (j24p is a differential consequence of another, simpler, equation 

3 A 

U^ e,,; - ^ ([/,,„ + Uui) = 0, (25) 

P P 

which is obtained from (|24p by differentiation and symmetrization with respect to the spatial variables. For a tidal 
action with metric (|23p. the boundary conditions have the form 

iXQN' + 2fiU.,kN'' = -X^^N' - f^^^N\ (26) 

Thus, in order to find the variations of the length of a rod due to variations of the gravitational potential, wc should 
solve homogeneous elastodynamic equations ()25|) with inhomogeneous boundary conditions p6|l . 

Under the above conditions for the gravitational field and for the elastic medium under investigation, there cannot 
exist a distinguished direction for the true deformations of the medium described by the tensor U ik ■ Therefore, we 
take the tensor of true deformations Uik in the form Uik = SikQ- 

For a one-dimensional rod of length Lq situated along axis x, the vibration equations and the boundary conditions 
are expressed as 

Ui = —XUi xx H — —Ui,xx, (27) 
P P 

C/i,, (0, t) = C/i,, (Lo, t) = (28) 
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In this case, the time variations in the variation of the rod length, Ai = U^{Lq, t) - U^{Q, t), are described by 



y£o I SLq 



^2^2 ^L(2m-1)2 ' 2m+l 



7r{2m—l)at 



7T{2m~l)at 



J (fCO 



7r(2m — l)ai 
' ~o 

7r(2m— l)at 



dt 



(29) 



dt]} 



X cos 



-7r(2m— l)x 



When 6(p{t) — Acos(figi), we have 
where r22m-i are the natural vibration frequencies of the rod. For the fundamental harmonic (m = 1), we have 



(30) 



(31) 



The length of the rod L(t) = Lq + AL(t) obtained under the above-mentioned assumptions can be calculated by the 
formula 



L (t) = Lo 



c2 7r2 c2 nl-Q^ 



3 J 



(32) 



In experiments on the verification of the LPI principle, one deals with a nonresonance case (i7i ^ ^g)', hence, formula 
(|32)) contains a very small coefficient of order fl'^/ill at a sufficiently small term S(p(t)/c'^. Therefore, in a practically 
important case for the experiment, we can neglect the extremely small terms and rewrite formula (|32p as 



L{t)= L 



fit) 



(33) 



Substituting these values into the formula flmo — t^ttic/ L(t) and neglecting the terms quadratic in (p(t)/c'^, we obtain 

'Pity 



(34) 



Using formula p6p . we obtain the following expression for the oscillation frequency of a laser without an absorbing 
cell in the resonator: 



1 + ^ 



(35) 



i.e., the phenomcnological parameter ^ in formula (jl7[) is equal to 1 



IV. A NULL GRAVITATIONAL REDSHIFT EXPERIMENT ON A TWO-RESONATOR LASER 

SYSTEM 

To resolve the above-mentioned alternative experimentally, one can use a stationary double laser system one 
of whose resonators contains a cell with a nonlinearly absorbing gas, which actually turns this resonator into a laser 
stabilized by a nonlinear absorption resonance. Figure 1 represents a possible optical scheme of the laser system (for 
simplicity, the electronic equipment necessary for tuning to the absorption peak of the cell is not shown in this figure) . 
The basic optical elements of this scheme (the cavity end mirrors 1, the gas-discharge tubes 2, the absorbing cell 3, and 
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the partially transmitting mirrors 4) are rigidly fixed on the common base 5. The mirror 6 and the semitransparent 
dielectric plate 7 provide the mixing of light radiated from two resonators of the laser system, and the photodetector 
8 serves for detecting the beat note. In a linear approximation in (? ■, one can obtain the following expression from 
(fT5)) for the oscillation frequencies (jJi.2 of the first and second resonators: 



^1,2 — ^0 



1 



1 



<5'l.2 



(36) 



where loq 



is the frequency of optical radiation, Si = 1 + pa - qa ^ ^ for the resonator without 
absorbing cell, and ^2 = 1 + {pa - Pb) - {<la - lb) ^1 for the resonator with the absorbing cell. 

A null gravitational redshift experiment by a double laser system consists in the following. Due to lunisolar tides, 
Newtons gravitational potential experiences a periodic time modulation, '^{t) = ^po + Sip{t), where 6ip{t) has a 
frequency of fig ^ 10"^ Hz and an amplitude of Sipo ~ 2.88 x 10^ m^/s^. Since the settling time of steady oscillations 



in the laser system is At <C27r/r2g {At is on the order of Af2^^ , where AQ^ 
(|36| to calculate the instantaneous oscillation frequency of the laser. 



10^-lO^rad/s), we can apply formula 
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— 7? 


'6 







Figure 1: Double laser system. 
According to (j36p . the difference frequency of two resonators contains a periodically varying term 



Auj 



(t) = ^„(i-0(^) 



(37) 



The value of this additional term essentially depends on the phenomenological parameter ^. For example, when ^ = 
0, the amplitude of the variations of the difference frequency Aa;( t ) for the optical radiation ( ~ 10^^ rad/s) is 
about 320 rad/s, which is much greater than the linewidth due to the spontaneous radiation of the atoms of the active 
medium in the case of gas lasers. Technical fluctuations of the difference frequency can be minimized by a single-block 
configuration of the double laser system and by placing it into a thermostabilized vacuum chamber shielded from 
external electric and magnetic fields. Moreover, since the detected signal is periodic, it can be accumulated during 
several months with the use of efficient methods for extracting a weak low- frequency signal from large noise [Toj . The 
potential accuracy of determining A/3 = — Pb\ (see formula ([3]) in the experiment is limited by natural fluctuations 
of the oscillation frequency of the laser that are associated with the spontaneous radiation of the atoms of the active 
medium and lead to a finite linewidth of the laser radiation. The latter parameter can be evaluated by the well-known 
Schawlow-Townes formula 

A.0 = ^ {A^nf , (38) 

where P is the power of laser radiation. For P \ mW and A^ji « 10^ rad/s, we obtain the following estimate for 
the least possible value of A/3: 



Af3 > l^-^j — (Anj^y « 10-'' (39) 
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which is four orders of magnitude higher than the accuracy achieved in the experiments of QiHl- 

Figure 2 represents another possible scheme of the detector that is based on a two-resonator laser system in which 
a common active medium is used for generating oscillations in two spatially inequivalent resonators. The first (signal) 
resonator is enclosed between mirrors R and Ri , and the second (reference) resonator is placed between mirrors R 
and R2 ■ The following optical elements are common to both resonators: the active medium AM (a gas-discharge 
He-Ne tube without Brewster windows), quarter- wave plates Li and L2 whose fast axes are mutually perpendicular 
and make an angle of 45° with the plane of the figure, and a WoUaston polarization prism P. The WoUaston prism is 
oriented so that the extraordinary ray with TM polarization (the electric field vector lies in the plane of the figure) 
is directed to the autocoUimation mirror Ri and the ordinary ray with TE polarization (the electric field vector is 
perpendicular to the plane of the figure) is incident to the autocoUimation mirror R2 ■ 



f-H-H 




Figure 2: Two-resonator laser system. 



First, we consider the operation of the two-resonator laser system in the absence of quarter- wave plates Li and 
L2- In this case, linearly polarized optical radiation with TE and TM polarizations is generated in the reference and 
signal resonators, respectively. The diaphragms Di and D2 serve for separating the fundamental transverse TEMgo 
modes in each resonator. In addition, by displacing one of the diaphragms perpendicular to the initial position of the 
optical axis of the corresponding resonator, one can significantly change the spatial overlap of the generated modes 
with TE and TM polarizations in the active medium; this sharply reduces the competition and coupling between these 
modes [2^. Radiation from the reference resonator emitted through the partially transmitting output mirror R2 is 
mixed, with the use of the beam splitter BS (see Fig. 2), with the radiation of the laser stabilized over the absorbing 
cell SL. The beat note detected by the photodetector PD2 is fed to the electronic unit of the automatic phaselocked 
frequency control device (PLFCD), which locks the oscillation frequency of the reference resonator to the oscillation 
frequency of the stabilized laser SL by means of a piezoelectric cell attached to the mirror R2 ■ Radiation from the 
signal and reference resonators transmitted through the common output mirror R passes through the linear polarizer 
Pq whose transmission axis makes an angle of 45° with the plane of Fig. 2 and forms an interference pattern, which 
is detected by the photodetector PDi . The signal from the photodetector PDl is fed to the signal-processing unit 
SP and to the electronic unit of the automatic frequency-control device (AFCD), which controls the operation of the 
piezoelectric cell attached to the mirror Ri. The AFCD enables one to fix the initial frequency difference between 
the signal and reference resonators or provides the operation within the frequency locking range when information on 
Auj{t) is transformed into the variation of the phase difference between the signal and reference resonators ^20]. The 
AFCD can also be used for partially compensating the noise due to technical fluctuations of the difference frequency 
of two generated modes with TE and TM polarizations. 

The minimal linewidth of the difference frequency of two laser modes is determined by the quantum noise level 
due to the spontaneous radiation of the atoms of the active medium. In 21| , Scully theoretically predicted that the 



linewidth of the difference frequency of two modes with different upper excited atomic levels and a common lower level 
can be substantially narrowed down (virtually to zero) by producing an active (compulsory) correlation of populations 
of the upper levels. In [2^.[23|. this idea was realized experimentally and was used to suppress quantum phase noise 
at the difference frequency of the Zeeman cr+ and (t~ components of optical radiation generated in a linear two-mirror 
single-mode He-Ne laser due to the splitting of the eigenmode of the resonator under a dc magnetic field Bq applied 
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along the active medium. The correlation of spontaneous radiation (the absence of the phase diffusion at the difference 
frequency) is achieved under the application of a transverse radio-frequency magnetic field B^foi a certain amplitude 
that couples the split upper laser levels due to a magnetic dipole transition between them. In the scheme shown in 
Fig. 2, such an effect can be realized by using quarter-wave plates LI and L2 under the application of a longitudinal 
dc magnetic field Bq and a transverse ac magnetic field Erf to the active medium AM. Applying the Jones matrix 
formalism to the investigation of the polarization state of eigenmodes in the two-resonator laser system, we can easily 
show that, as before, optical radiation with TE and TM polarizations is generated to the left of the plate Li and 
to the right of the plate L2 in the reference and signal resonators, respectively, which correspond to clockwise and 
counterclockwise polarized waves in the space between the plates. In this operation mode of the detector, the time 
needed to extract a useful signal can be significantly reduced because of the absence of a phase drift at the difference 
frequency, which is associated with the spontaneous radiation of the atoms of the active medium. 

V. CONCLUSIONS 

Comparing formulas (1) and (11), we can easily see that the quantity A/3 =\Pa — Pb\, which characterizes the 
deviation from the LPI principle in a null gravitational redshift experiment, restricts the possible values of the 
phenomenological parameter to |1 — ■CI ^ ^f^- The experiments of '?|,'^ show that A/3 is no greater than 10~^; 
therefore, the parameter . must be close to unity. In the case of an experiment carried out according to the scheme 
proposed in this paper, one can obtain more accurate values of A/3 (to within lO^^or less). This would allow one to 
verify the LPI principle for clocks of different physical natures more accurately than was done in At the same 

time, one may experimentally verify the correctness of the relativistic generalization of the classical elasticity theory 
considered in this paper. 
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